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Abstract. - We predict theoretically that long- wavelength surface charge modulations uni- 
versally reduce the pressure between the charged surfaces with counterions compared with the 
case of uniformly charged surfaces with the same average surface charge density. The physical 
origin of this effect is the fact that surface charge modulations always lead to enhanced counte- 
rion localization near the surfaces, and hence, fewer charges at the midplane. We confirm the 
last prediction with Monte Carlo simulations. 



Introduction- The interaction of charged interfaces is a basic problem in surface science, 
and is relevant to e.g., colloid stabilization, membrane adhesion, microemulsion formation 
While many theories have focused on the simple case of homogeneously charged surfaces, real 
interfaces are characterized by discrete surface charge distributions. Theoretical analysis of 
the problem predicts two important properties of these systems: the spatial dependence of 
the counterion/salt density and the pressure between the interfaces fl|3|,fl|5|,|6|,[7||^,p|, [l0||ri] ,|l^, 
The main question is whether one can provide universal predictions of the effect 
of surface charge discreteness (modulations) on the counterion spatial distribution and on the 
interaction law between the interfaces, independently of the details of the surface distribution 
of discrete charges (modulations). Our work shows that indeed, there are quite general effects 
induced by surface charge discreteness. 

The principal conclusion of previous studies based on the solution of linear or quasi- 
linear Debye-Huckel equation with non-uniform surface charge distributions with cither excess 
salt H,^,^, or with counterions only jl^], is the following: The discreteness (or modulations) 
of surface charge always leads to an enhanced repulsive interaction between the surfaces com- 
pared with the case of uniformly charged surfaces with the same average surface charge density, 
provided the charges (or charge modulations) on the two surfaces are in-phase. If the charges 
(charge modulations) on the two surfaces are out-of-phase, the repulsion can be either en- 
hanced or reduced compared with the uniformly charged surface case. These studies [p|p|,p^ 
calculated the effects of charge modulations on the counterion spatial distribution and on the 
electrostatic potential to first order in the amplitude of modulations. The only relevant terms 
in the pressure are therefore due to the square of the non-zero components of the electric field, 
E, at the midplane ^E'^-\- Ey — El (see below). This contribution is always positive for the 
in-phase case, and leads to the enhanced repulsion in the in-phase case. 
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In this paper, we show that a consistent calculation of the pressure, including both the 
osmotic terms proportional to the counterion charge density at the midplane, and the electric 
field terms calculated to second order in the amplitude of modulations, results in a universal 
reduction of the pressure by long-wavelength surface charge modulations. The notion of a 
universal reduction means that the effect does not depend on a particular form, amplitude, and 
phase shift between surface charge modulations. The physical origin of the predicted reduction 
of the pressure is due to the osmotic effect and arises because surface charge modulations 
always enhance counterion localization in the vicinity of surfaces, and hence, always reduce the 
midplane counterion density. We confirm this prediction for the charge density with Monte- 
Carlo (MC) simulations. We emphasize that the predicted effects are due to the correlations 
between the counterions and the inhomogeneities of the surface charge distribution. 



Two surfaces with charge modulations- In this section we consider two, non-uniformly 
charged planar surfaces with fixed surface charge densities —eai{p) and —ea2{p) [where p = 
(x, ?/)], respectively. The surfaces whose average charges are equal and are taken to be negative, 
are located at 2; = — | and 2 = f , respectively. The surface electrostatic repulsion is screened 
by the positively charged counterions of valence Z; the counterions are located only in the 
water between the surfaces, and the dielectric constant, e, (e ~ 80 for water), is assumed 
to be homogeneous throughout the system. The charge densities cri(/5) and <J2{p) obey the 
condition: / a-i{p) dp/Ao = ao, where Aq is the surface area, and ao = Z N / 2Ao is the average 
number of surface charge per unit area, with N being the total number of counterions in the 
system. 

The system is fully specified within mean-field theory, by solving the differential equation 
for the potential |p^, ip(p,z), with the boundary conditions implied by the delta function, 
spatially dependent surface charge distributions ai (p) and a2 {p) ■ 



V'^4>{p,z) = -AttIbZ 



Zno e-^ - 5{z + ^)a,{p) - 5(z - ^)a2{p) 



\A < \, (1) 



and {p, 2) = A'k£bZ[S{z + |)cti(p) -I- S{z — |)ct2(p)], for |2| > ^, where there are no coun- 
terions ||l6|. Here 4>{p,z) — eZLp{p, z) /ksT is the reduced electrostatic potential, — 7^7^ 
is the Bjerrum length, and 5{z) is the Dirac delta function. In what follows, we con- 
sider a general, modulated charge distribution described by a Fourier decomposition [nJ: 



(JiiP) = cro 1 + EQ^o^('5)e"^ '^ and (72 (p) = ao 1 + Eq^^o '^(Q) « 



tQ-ip+b) 



where & is a 



phase shift of the modulations on the two surfaces. The coefficient e{Q) is the amplitude of 
the corresponding Q mode; for example, on a square lattice, Q — 2Trrn/a, and m — {mx,my), 
where rux and my are integers, and a is a lattice constant. For a periodic lattice of sur- 
face charges, one must keep the infinite sum over reciprocal lattice vectors, Q, and set 
£{Q) = 1. In analogy with the analysis performed earlier for the case of a single, isolated 
surface we expand the potential, 4'{p, 2), in powers of the charge modulation amplitude, 
e{Q): (j3{p,z) K, (j)o (2) -)- 0i (/?, 2) -1-02 {p^z) + where 00 {z) — 21ncosfco2 is the solution of 
the PB equation with uniform surface charge densities, —etro. The parameter fco is deter- 
mined by the charge conservation condition jl^: fcoAotan(^Y^) = 1, where Ao = l/^irlBZao- 
Our objective is to obtain perturbatively the corrections 0i (p, 2) [linear in e{Q)] and 02 (p, 2) 
[quadratic in e{Q)]. Substituting the expansion of (j){p,z) into Eq. (|l|), and collecting terms 
linear in e{Q), we find that (pi{Q, 2) [the two-dimensional Fourier transform of 0i {p, z)] obeys 
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the following equation for \z\ < 

2eiQ) 



2 fcp 



dz-^ cos^ koz 



An 



8{z + ^) + 5{z-^)e^^-^' 



, Q ^ 0. (2) 



From here on, all the lengths in the problem will be rescaled in the units of Aq: Q = QAq, 
z = z/Xq, p = p/Xq, ko = fcoAo, h = h/Xo, etc. The solution of this equation, 0i( Q, z), is 
given in the Appendix. 

The equation for the contribution to the potential that is second order in the surface 
charge modulation amplitude, (j)2{p,z), is obtained in a similar manner; collecting the terms 

quadratic in e(Q) we find for \S\ < |: f-V^ + -J^l Mp,S) = T^Jfe This 



equation is solved, using the Green function for the homogeneous equation, |18|. In what 
follows, for the sake of simplicity, we consider two opposite cases: (i) in-phase, i.e., maximum 
interaction, b = (0,0), and (ii) out-of-phase, i.e., minimum interaction, b = (|, |-). The case 
of arbitrary b may be solved analytically, as well. We will consider surface charge modulations 
described by a single symmetric Q mode: Q = {(±^,0), (0, ±^)}. 

The average of 02 in the x — y plane is given in the Appendix in both cases. Our main 
prediction is the enhancement of the counterion density at each surface, averaged over the x—y 
plane, {n{p, z = ±-|)) _.. The physical importance of this result is that by charge conservation 
it leads to the reduction of the midplane average density compared to the case of uniformly 
charged surfaces in both the in-phase, b — (0,0), and out-of-phase, b = (§,§), cases. This 
conclusion comes from the analysis of {n{p, 5)) - which takes the following form to quadratic 

order in e(Q): (n(p,5))- = no(z)[l + \ (^^(p, z))^ - (02(/9,z))-], where n^iz) = hq/ cos^ koZ 

is the counterion density between uniformly charged surfaces fcg — 27r^BZ^no. We 

emphasize that the first-order term, (f)i{p,z), vanishes after averaging with respect to p, and 
that only quadratic terms in e(Q) contribute to the average density; this observation is crucial 
for the correct calculation of the osmotic pressure, as we show below. The enhancement of 
the counterion density in the vicinity of the surfaces [and the consequent reduction of the 
midplane density] is due to the correlations between the counterions and the inhomogeneities 
of the surface charge distribution. In particular, in the asymptotic limit of small inter-surface 
separation, h <^ I, Qh -C 1 the normalized average midplane counterion density has the 

form: (i) in the in-phase case, ^"/.^ff)^'^ « 1 - |I]g (2+^q)=' ^ + ^i^^)y (") ^^'^ out-of- 
phase case, ^"^^(0)^*^ ~ ^ ~ ^ J2q + 0(^"^)- We stress that the strongest reduction 
of the midplane counterion density [and hence, the strongest enhancement of the contact 
density] occurs for the long-wavelength modulations. This is in a qualitative agreement with 
the predictions obtained for a single, isolated surface [Q. In the opposite limit of large 
separations, /i 3> 1, Qh 3> 1, the midplane density is reduced with an asymptotic form that 
is identical in both limits, (i) and (ii): ^^g^ « 1 " Eg wKtW^ + 

To verify the predictions of the theory, we have performed MC simulations of the counte- 
rion density profiles between non-uniformly charged surfaces. The non-uniform surface charge 
densities were chosen to have the form: ai{p)/ao = 1 -f eoX^g e*^'^, and a2{p)/cro = 1 + 
So X]q ■ip+f>)^ respectively, where Q takes on the four symmetric values: {(±^, 0), (0, ±^)}. 
The counterions interact via an unscreened Coulomb potential, and an Ewald 2D algo- 
rithm [p^ was used to sum the Coulomb interactions; each counterion interacts with the 
surfaces via the potential, eZips generated by the exact solution of the Poisson equation with 
the chosen surface charge densities: ^ = -eo Eg ''^(6"®'^+''/^' + e*^ ■be-2l^-''/2|)/Q. 
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Fig. 1 - Normalized counterion density in the in-phase case, n{z) = 27r^sZ^Ao {n{p,z))^ from the 
MC simulations (symbols) and the corresponding predictions of the analytic, perturbation theory 
(curves). The parameters used are: g = — 0.5, eo = 1, Q = 1 (diamonds and dashed line). 

The counterion profile between uniformly charged surfaces with the same average surface densities 
of charge (stars and solid line) is shown for comparison. The simulations were performed with 100 
particles. Inset: The pressure difference diagram, p — po = 0, as a function of the rescaled wavelength 
of modulations, Q, and the inter-surface separation, h, in the in-phase case, where po is the pressure 
between uniformly charged surfaces. We stress that already at /i > 1.25, charge modulations reduce 
the pressure for any values of Q. 
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Using the simulations we compare the counterion density averaged in x — y plane. We plot 
the normahzed dimensionless density, n{z) = 2ttIbZ'^^o {^{P:^))^ vs. the rescaled length, 
z = z/Xo- The analytic expressions for {n{p, z))~ are given above (see also the Appendix) for 
both the in-phase and the out-of-phase cases. 

The dimensionless parameter which determines the strength of the Coulomb interaction in 
the simulations is g = ^^'^ ; other two important parameters are the amplitude of the surface 
charge modulation, eo and the lengthscale, Q, of the surface charge modulations. In Fig. 1 we 
plot the simulated, average density, n{z), in the in-phase case, b = (0, 0) for values of g — 0.5, 
Eq = 1, and Q = 1 (diamonds), as well as the simulated counterion profile between uniformly 
charged surfaces with the same average surface charge density (stars). Additional simulations 
showed that the average density is higher for the smaller values of Q, similar to the conclusion 
obtained earlier in the case of a single, isolated surface We also note that the quantitative 
difference between average density profiles predicted by theory and confirmed by simulations 
in the in-phase and the out-of-phase cases is very small. The contact (midplane) density is 
always slightly higher (smaller) in the in-phase case compared with the out-of-phase case, since 
in the out-of-phase case the broken planar symmetry leads to the existence of an oscillating 
z— component of the electric field at the midplane (absent in the in-phase case), that attracts 
more counterions to the midplane due to the enhanced correlations. We stress again that 
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the perturbation theory prediction of the enhancement of the contact counterion density is 
universal in the sense that the enhancement occurs for all values of the parameters, i.e., Q, h, 
and £, in both the in-phase and the out-of-phase cases. We emphasize that ion fluctuation and 
correlation effects neglected in the mean- field approach we adopt, become important at high 
values of 5 ^ 1, i.e., in the strong coupling regime jT^ . Finally we note, that qualitatively 
similar conclusions with respect to the enhancement of the contact counterion density have 
recently been obtained by numerical solutions of the integral equation theory p^ ; however 
some of their effects are related to the fact that both the counterions and the surface charges 
have an excluded volume in [ pl| . 

Pressure between charge modulated surfaces- In this section, we show why surface charge 
modulations in the in-phase case (under certain conditions), and in the out-of-phase-case 
(for all values of the parameters) always reduce the repulsive pressure between surfaces with 
counterions compared with the case of uniformly charged surfaces. 

The free energy, F, of the system has the form = J n(r )[ln(n(r )wo) — 1] -f 

1/ '/>(^)'T'tot(?') rf?', where the total system charge density: ntot{r) — n{f) — ^[ai{p)6{z + ^) + 
a2{p)S{z — ^)], Wo is an ion-size volume, and r = {p,z). Keeping terms up to second order 
in the amplitude of surface charge modulations, s{Q), one notes that F can be expressed 
only in terms of the x — y average of the potentials as obtained above, {4>2{r ))^ , i't'ii'^)) p 
and {o'i{p)4>i{r)) after straightforward but lengthy integrations, one obtains the following 
dimcnsionless free energy, / = k^TAoia-o/z) ' P^'^ counterion /_(- and /_ in the in-phase and in 
the out-of-phase case, respectively: 



where fo is the free energy of uniformly charged surfaces with counterions . The normal- 
ized, dimcnsionless osmotic pressure between the surfaces, p = —df/dh, is the main quantity 
of interest measured in experiments. Our principle observations follow intuitively from the 
asymptotic analysis of the pressure: (i) in the in-phase case, in the limit of small separations, 
h <^1, Qh <^ 1, one has. 







where po — is the normalized, dimcnsionless pressure [pO| between uniformly charged sur- 
faces jl^. In the out-of-phase case, when /i <C 1, Qh ^ 1 we obtain, p — po = — X] q + 
0{h). Therefore in this asymptotic limit of small inter-surface separations, in the in-phase case 
the pressure between the surfaces is reduced for long- wavelength surface charge modulations; 
in the out-of-phase case, the pressure is reduced for any wavelength. This effect [quadratic in 
e(Q)] is related to the reduction of the midplane counterion density induced by surface charge 
modulations upon decreasing of Q. In the limit of large inter-surface separations, /i 3> 1, 
Qh 3> 1, in both cases, (i) and (ii), the pressure is always reduced: 





Q 



We note that this asymptotic form is similar to the universal, Casimir-type attractive fluctua- 
tion pressure, —\/h^ (see e.g., [^|2^) between the surfaces with boundary conditions that 



6 



EUROPHYSICS LETTERS 



suppress or modify fluctuations of the medium; the amphtude of this pressure is non- universal, 
of course, since the spectrum of modes is restricted by the fixed surface charge modulations 
in the present case. 

It is interesting to note, that corresponding asymptotic results for the pressure pO| ] between 
uniformly charged surfaces po [ po — j^, if h ^ 1, and po — f^, if h ^ 1] show that 
the modulation correction, p — po, can never overcome the overall repulsion in these limiting 
cases. For intermediate values of h, however, the total pressure, p, can change sign and become 
attractive, even in the in-phase case, if we extrapolate our perturbative results to values of 
e(Q), which are not necessarily small. We also emphasize again that we use a mean- field 
approach, and neglect the fluctuation contribution to the pressure | pT| that further reduces 
the repulsive interaction. 

The generic behavior of the pressure difference, p — po, in the in-phase case for arbitrary 
h and Q is shown in the inset of Fig. 1. There is a broad range of parameters where the 
repulsion between two surfaces is reduced by the effect of surface charge modulations. It is 
interesting to note, that already at /i > 1.25, charge modulations reduce the pressure for any 
value of Q. We have checked that a qualitatively similar reduction of the repulsive pressure 
occurs in the case of a square lattice of discrete charges (infinite sum with respect to all Q 
modes). 

Our prediction of the reduction of pressure in the in-phase case differs from the con- 
clusions of Rcfs. @,||l|. The differ ence between those works and our predictions can be 
understood using an alternative definition of the pressure via the stress tensor p = 

[n+ J + (j)y — (^1) _.]i=o , where n = 2tt£bZ'^^o (n)^ is the normahzed, dimensionless den- 
sity |2^; this gives results identical to the ones presented above. In the in-phase case, (jj^ = 0, 
as implied by the symmetry, and the electrostatic contribution to the pressure is purely re- 
pulsive, arising from terms in the electric field that are first order in the amplitude of sur- 
face charge modulations, e{Q) [the pressure is, of course, always quadratic in s{Q), since 
p (fi^ + (j)y]. This is the only contribution that was taken into account in Refs. P, p|p^, 
where only first order corrections to the counterion density were calculated. As far as the 
osmotic term is concerned, the first order corrections vanish when the average over the x — y 
plane is performed, and hence, it is necessary to calculate the counterion density to second 
order in the surface charge modulation amplitude, e{Q). It is this osmotic contribution, 
calculated here, that is the origin of the universal reduction of the pressure, since n at the 
midplane is always reduced by surface charge modulations, as we have shown; this effect is 
quadratic in e(Q), and beyond the linear approach of Refs. ||,^,|l5|. It is the competition 
between the osmotic and electrostatic contribution that leads to the diagram represented in 
the inset of Fig. 1. In the out-of-phase case, where 4>s ^ 0, the electrostatic contribution to p 
is reduced compared with the in-phase case, and further reduction of the midplane density by 
charge modulations always leads to the reduction of the pressure, compared with the case of 
uniformly charged surfaces. 

In summary, we emphasize that our main result, the prediction of a universal reduction of 
the osmotic pressure by long-wavelength, quenched surface charge modulations, is naturally 
related to effects of fluctuating surface charges on the pressure. The unifying theme is that 
systems with thermally fluctuating surface charge density are qualitatively similar to systems 
with quenched, modulated surface charge density: in both situations, the counterions are more 
localized near the surfaces compared with the case of uniformly charged surfaces, and in both 
cases this leads to a reduction of the pressure. 



* * * 
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Appendix- In this appendix, we summarize the results for perturbative corrections, 4>i{p, z), 
and {(f>2{p, z)) ^, to the potential, 4>o{z), induced by surface charge modulations. 

In particular, one finds: z) =^ Avi(Q,S) + Bv2(Q,S), where ^ = M^l "i+g2e' " ^ 

g = «i = exp(- fi)(fc2 + i)/Q2^ ^ exp(fi)(fcg + 1 + 2Q + 2Q^)/Q^ 

and vi{Q,i) = e^''{l + {ko/Q) tanfcoi), f2(S,5) = e-2^(l - (fco/Q) tanfcoi). 

The average potential, {4>2{p,z))^ is also obtained in a simple, closed form in terms of 

elementary functions: {4'2{p, z)) ^ ^^cii^^^|=pat j^^ ' where r+(Q) and r_(Q) corresponds 

to the in-phase case and the out-of-phase case, respectively. Here we use the notations: 
r±(Q) = U2{i)[^g±{^) - ^±(Y) + -^±(^0^)] - u^{~z)g±{h~z), u^{i) = tanfcoz, U2{~z) = 
l+fcoztanfco5,r; = ^ + ^,e±(a;) = f|h||^ [tanx-x(l-tan2 x)]H-fi|f| [1 a;tana;]± 

1 [ ^ + tanxl , and T±{x) = [^"'^h^e^Tl] tan^. sinh2ex tan. ^ 1^ _ co^ 
4 L cos^ X J ' =n y / g ^2 4y cos^ a: 4 cos^ x 8 ' 

= Q/fco- 
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